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We investigate the BCS-BEC crossover in three dimensional degenerate Fermi gases in the presence 
of spin-orbit couphng (SOC) and Zeeman field. We show that the superfluid order parameter 
destroyed by a large Zeeman field can be restored by the SOC. With increasing strengths of the 
Zeeman field, there is a series of topological quantum phase transitions from a non-topological 
superfluid state with fully gapped fermionic spectrum to a topological superfluid state with four 
topologically protected Fermi points {i.e., nodes in the quasiparticle excitation gap) and then to 
a second topological superfluid state with only two topologically protected Fermi points. The 
quasiparticle excitations near the Fermi points realize the long-sought low-temperature analog of 
Weyl fermions of particle physics. We show that the topological phase transitions can be probed 
using the experimentally realized momentum resolved photoemission spectroscopy. 
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Ultra-cold Fermi gases with tunable atom interaction 
through Feshbach resonance [l| have garnered tremen- 
dous attention recently for their potential use as 
an ideal platform for emulating many important phys- 
ical phenomena. Interesting physics, including the 
crossover from Bardeen-Cooper-Schrieffer (BCS) super- 
fluid to Bose-Einstein condensate (BEG) of molecules ^3], 
universal properties in the unitary limit vortices [5[, 
etc., have been observed in experiments. In addition to 
single species Fermi atoms with equal population in two 
spin states, Fermi gases with population and mass im- 
balance have also been intensively studied (Jll]. Here 
the population imbalance between two pseudo-spin states 
serves as an effective Zeeman field, which have stimulated 
enormous experimental efforts on searching for the Fulde- 
Ferrel-Larkin-Ovchinnikov state in Fermi gases 0. 

The pseudospin of atoms can couple with not only the 
effective Zeeman field, but also with the orbital degrees 
of freedom of atoms {i.e., spin-orbit coupling (SOG)). 
The recent experimental realization of SOG for ultra- 
cold atoms 10] opens a completely new avenue for in- 
vestigating the SOG physics using degenerate cold Fermi 
gases. In this context, it is natural to investigate the 
BGS-BEG crossover physics 11, [l2] in the presence of 
SOG, which is an important and interesting problem by 
itself and has not been investigated previously. Further- 
more, such crossover physics is also important because 
the s-wave superfluid, together with SOG and Zeeman 
field, may yield intriguing chiral p-wave physics [13] such 
as Ma iora na fermions with non-Abelian statistical prop- 
erties |14l |. However, such chiral p-wave physics may be 
observable only in the crossover region where the super- 
fluid order parameter is large and thus robust against 
flnite temperature effects. 

In this Letter, we investigate the BGS-BEG crossover 
and the topological properties of a three-dimensional 
(3D) uniform s-wave superfluid in the presence of both 



Zeeman held and Rashba-type of SOG. Under the mean 
held approximation, we derive the superfluid gap and 
atom density equations and solve them self-consistently 
in the BGS-BEG crossover region. Our main results are: 

(I) It is well known that the superfluid can be destroyed 
by the Zeeman field beyond a critical value for a given 
s-wave interaction strength ^5i]. We show that a finite 
SOG strength can restore the superfluid pair potential 
back to the system even when the Zeeman field is well 
above the critical value. 

(II) At zero temperature, the nonzero superfiuid pair 
potential in the presence of SOG supports, with increas- 
ing strength of the Zeeman field, a series of 3D topological 
quantum phase transitions from a non-topological su- 
perfiuid state with fully gapped fermionic excitations to 
a topological superfiuid state with four protected Fermi 
points {i.e., gap nodes) and then to a second topologi- 
cal superfiuid state with only two protected Fermi points 
(see Fig. 3). Such 3D topological superfiuids are gapless 
without Majorana fermions because of the existence of 
Fermi points, as opposed to the 2D fully gapp ed topo- 
logical superfiuids with Majorana fermions 14 1. 

(HI) The superfluid phases separated by the topologi- 
cal quantum critical points are indistinguishable in terms 
of the superfluid pair potential, which remains continuous 
across both transitions. On the other hand, the super- 
fluid phases are distinguishable in terms of the ground 
states which have very different excitation spectra. The 
non-topological superfluid state is a usual fully gapped 
s-wave superfluid, while the topological superfluid states 
are gapless even if the superfluid pair potentials are still 
strictly s-wave! These peculiar gapless topological super- 
fluids are new phases of matter and here we show how 
to identify such phases and the corresponding quantum 
critical points using the experimentally already realized 



momentum- resolved photoemission spectroscopy |17l.ll 



Mean-field theory: The system we consider is a 3D uni- 
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form s-wave fermionic superfluid with the atom density 
n and which is subject to Rashba SOC in the xy plane 
and a Zeeman field along the z direction. The dynamics 
of the Fermi gas can be described by the Hamiltonian 
H = Hq + Hint, where the single particle Hamiltonian 
Ho = Ek77' ^7 Kk/ + a (fcyCr^ - fcxCTy) +rcr2]^^, Cky, 
7 -I are the pseudo-spin of the atoms, £,-k — e-k — 
Ek = h^k"^ /2m is the free particle energy, ^ is the chemi- 
cal potential, a is the Rashba SOC strength, / is the 2x2 
unit matrix, ai is the Pauli matrix, T is the strength of 
the Zeeman field, and Ck7 is the atom annihilation opera- 
tor. i?int — g X^k '-kt'^-ki'^-ki'^kt is the s-wavc scattering 
interaction with g — Awfi^as/m, and the scattering length 
fls can be tuned by the Feshbach resonance 

In the mean-field approximation, the 5-wave super- 
fiuid pair potential A — g (ck4,c_kt} ^nd -ffi„t — 
-A^/g + AX]k(ckiC-kt+ cLkt4i)- Under the Nambu 
spinor basis ^Pk = (ckti Ckj., cLj^^, — cLj^^)^, the Hamilto- 
nian is if i Ek *lc^^k*k - AV.9 + Ek ^k, where 
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preserves the particle-hole symmetry. The quasiparticle 
excitation energy 



r2 + |A|2±2So 
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is the eigenvalue of Afk, where A = ± corre- 
spond to the particle and hole branches, Eq = 

VrW+WT+^^^, and ki_ = ^kl + k 

a = r = 0, Eq. @ reduces to = Xy^^+\A^ in 
the standard BCS theory. 

Diagonalizing the total Hamiltonian using the Bogoli- 
ubov transformation and following the standard proce- 
dure we obtain the gap equation 
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where / J = tanh(/?i?+^/2)/4i?+„, /3 = l/fcsT, T 

is the temperature, and ks is the Boltzmann constant. 
Following the standard procedure pi] , the ultra-violet di- 
vergence at the large k in Eq. ([3]) has been regularized 
through subtracting the term 1 /2ek and is defined as 
the renormalized scattering length. The total number of 
atoms can be obtained using a similar method [2] 
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We self-consistently solve the gap equation ^ and 
the number equation ([4]) for different parameters 
(tt-ftTF, r, z/, T) for a fixed atom density n to determine 

A and /i. Here v = l/Kpa-s and Kp = (Stt 



'^n]^^^ is the 
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FIG. 1. (Color online) Plot of /i (a) and A (b) versus v 
l/Kpas for different parameters [aKF,!") at T = 0. 



Fermi vector for a non-interacting Fermi gas with the 
same density at F = a = 0. The energy unit is chosen as 
the Fermi energy Ep = h?Kp/2m. 

BCS-BEC crossover: Without SOC, it is well-known 
that the Zeeman field F can lift the spin degeneracy be- 
tween t and 4,, thus destroys the Cooper pairing when 
the Zeeman field is larger than the pairing interaction 
energy flst . This result can also be understood from 
Eq. ([5]), where the quasiparticle excitation gap Eg = 

= for a suitable chosen k when 



lAP-F 



mink Vl^ 

a = and |F| > A. With SOC, the single particle Hamil- 
tonian Hq has two bands and each band contains both 
spin 'I' and 4, components even with a large Zeeman field, 
leading to nonzero superfluid Cooper pairing between two 
fermions in the same band with opposite momenta. 

Such SOC induced nonzero superfluid pair potential 
can be clearly seen in Fig. [1] where we plot the change 
of and A with respect to v = l/KpUs for different 
parameters {aKp,T) at T = 0. In the BEC side, the 
SOC strength aKp and the Zeeman field F do not have 
significant influence on A because all fermion atoms form 
bound molecules. Hence we mainly focus on the BCS 
side. At a = 0, the superfluid pair potential is destroyed 
when F > A, as expected. In contrast, A and fj, are 
independent of F when F < A. This can be understood 
from the fact that, without SOC, Eq. ^ and (g]) are 
independent of F when F < A. Therefore there is a 
sudden jump of fi and A at F = A for a — 0, as clearly 
seen from Fig. [TJ We see A can be restored to non-zero 
values even for a large F when aKp is switched on. 

In Fig. [5^, we plot A with respect to aKp on the 
BCS side with i/ = —1 and T = 0. For other values of 

and T, the results are similar. When F < A, the pair 
potential A approaches the same point (determined by 
the interaction strength v) for different F when a — s- 0. 
A vanishes for F beyond a critical value at a = 0, but can 
be restored when a is nonzero. Therefore the superfluid 
order can still be observed even with a large Zeeman flcld. 
In Fig. [2b, we plot A versus F for different SOC. Without 
SOC, we observe a sudden jump of the pair potential at 
F = A, as expected. With non-zero SOC, A decreases 
smoothly with F. At large F, the numerical results can 
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be fitted with A ^ ^ witli tlie constant x depending 
on the SOC. 

Topological phase transition: The existence of nonzero 
superfluid pair potentials induced by the SOC does not 
automatically determine the topological properties of the 
superfluid. In a 3D uniform superfluid, the momentum 
is a good quantum number and the topological order of 
the superfluid can be classified by the number of Fermi 
points (i.e., the nodes such that _ = 0) in the 
quasiparticle excitation spectrum for fixed external pa- 
rameters (aKpjT,!^) and at T = 0. From Eq. 1^, we 
seeE+_E+^ = {^l-V^ + A^-a^klf+AA^a^kl, there- 
fore there is always a finite gap for fc_L > and all Fermi 
points must lie along kz axis, see Fig. [SK- At k± = 0, 
the quasiparticle excitation gap becomes 
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FIG. 3. (a) Fermi points along the kz axis, (b) Quantum 
phase transition in a 3D uniform superfiuid driven by the 
Zeeman field at T = 0, aKp — Ef, and v — —0.1. The solid 
line represents Eg = 0. The change of Eg along the red arrow 
is shown in Fig. [4^. 
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where p, = /j, — h k^/2m. 

We choose aKp = Ep, v — —0.1, T = 0, and study 
the topological phase transition driven by the varying 
r. In Fig. ISb, we plot -Eg = in the {k^.V) plane 
and count the number of Fermi points for a fixed F. 
There are three different phases: (A) a fully gapped non- 
topological superfluid state without Fermi points when 
F < A, (B) a superfluid state with four Fermi points 
for each F when A < F < F^ = and (C) a 

second superfluid state with two Fermi points for each F 
when F > Fc. Around each Fermi point, we have checked 
that the quasiparticle energy dispersion E^_ is linear, 
therefore the Fermi point behaves like a Dirac point, 
and is topologically protected with a topological charge 
Na = ±1. Here Na is defined by the 3D topological 



invariant Na — 
[l6j , where = ik^ 
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Mk is the Green function for the 
quasiparticle excitation, is a 3D surface around the 
isolated Fermi point (0, k^), fco is the energy, and "tr" 
stands for the trace over the relevant spin indices. The 
quasiparticle excitations near the Fermi points realize the 
long-sought low-temperature analog of Weyl fermions of 
particle physics. At each of the two critical points as 
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FIG. 2. (a) Plot of A with respect to aKp for different F 
(b) Plot of A with respect to F for different aKp. v — ~1 
T = 0. 



points with opposite topological charges merge and form 
a topologically trivial Fermi point with Na = 0, which 
then disappear with the Zeeman field. Note that topo- 
logically protected Fermi points with Na = ±1 in phases 
B and C are protected in the sense that they cannot dis- 
appear with F except at quantum critical points (where 
they pair-wise merge to become topologically trivial and 
disappear with the Zeeman field |16|). Such a series of 
two topological quantum phase transitions (two quan- 
tum critical points) is quite different from the situation 
encountered across a p-wave Feshbach resonance (i.e., 
driven by varying the interaction strength), where a 3D 
gapless px + ipy BCS superfluid becomes a gapped BEC 
at a single topological critical point at ^ = [3]. In 
addition, the superfluid pair potential is still finite in our 
system, while it vanishes in the px + ipy superfluid [l9l | 
when the quasiparticle excitation gap closes. 

Observation of topological phase transition: The pe- 
culiar phases B and C with strictly s-wave pair poten- 
tials, but topologically protected Fermi points, are new 
phases of matter. Here we show how to identify such 
phases and the corresponding topological quantum crit- 
ical points by the momentum resolved photoemission 
spectroscopy [13, Note that topological quantum 

critical points are not easily observable in usual experi- 
mental probes because of the absence of broken symme- 
tries and emergent order parameters at such transitions. 
In the photoemission experiments, E^ _ for a flxed (or 
a small range of) kz can be measured by analyzing the 
part of the time of flight image with the same kz- This 
measurement yields Eg as a function of kz and F, from 
which the phase diagram in Fig. [3}d can be mapped out 
through Eg {kz,T) = 0. In the following, we describe the 
change of Eg for kz ~ and increasing F (the horizontal 
red arrow in Fig. [3}d). The results are similar for kz ^ 0. 

For kz — 0, Eg does not close at the boundary between 
phases A and B, therefore the quasiparticle excitations 
are gapped non-topological when F < Fc (i.e., A and 
B phases merge together). Eg vanishes at Fc and then 
reopens when F > Fc, where the quasiparticle excitations 
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FIG. 4. (a) Plot of Eg, A, and fi, with respect to F at k = 0. 
v = —0.1, aKp = Ef, and T — 0. (b) Topological phase 
crossovers at finite temperature, where the red curve is deter- 
mined by Eq. (6) The arrow represents the phase crossovers at 
a fixed temperature. A and C represent regions with gapped 
non-topological and topological excitations, respectively. 



become topologically non-trivial. In Fig. S^, we plot Eg 
at k = with respect to F at T = 0. The corresponding 
A and /i are also plotted. When F sweeps through Fj, — 
0.56Ep, Eg first closes and then reopens, indicating a 
transition from non-topological quasiparticle excitations 
to gapless excitations and then to topological excitations. 
At Fc, we still observe a strong A = 0.52Ef, although Eg 
has a node. For other fixed kz ^ 0, the corresponding Fc 
can also be determined similarly through the vanishing of 
Eg. Therefore the phase diagram for the 3D topological 
phase transition in Fig. [3b can be mapped out. 

A realistic experiment only works at a finite tempera- 
ture, which induces thermal excitations and may destroy 
the superfluid. The percentage of the thermal component 
in the gas depends strongly on the ratio E'g/fcsT. Note 
that at finite temperature, there is no definite boundary 
between different phases, i.e., there are only finite tem- 
perature phase crossovers. Nevertheless, we still plot the 
phase lines in Fig. |4}d at which 



Eg {a,i/,r,T) = ksT 



(6) 



for the illustration of the topological phase crossover re- 
gion. Note that Eg also depends on T, and Eq. © 
should be solved self-consistently with Eq. ([3]) and (|4]). 
Therefore the relationship between A and F here is quite 
different from the zero temperature case in Fig. |4^. In 
experiments, the fixed temperature corresponds to a hor- 
izontal line (with the arrow) in the phase diagram Fig. 
|4}d. Between phases A and C, there is a quantum critical 
region where finite temperature generates many quasi- 
particle excitations. 

Experimental parameters: In experiments, the Rashba 
SOC and Zeeman field can be generated through the 
coupling between atoms and laser fields [13, . We 
consider a set of realistic parameters for atom with 
n = 5xl0^^ cm~^, corresponding to the Fermi energy 
Ef = hx 3.5 KHz. A large SOC strength aKp ^ Ep 
is also attainable [lol . [20| . The effective Zeeman field F 
from to Ep can be obtained by tuning the intensity 



of laser beams 13, 2^ • By tuning the s-wave scattering 
length using the Feshbach resonance we can vary A 
from ^ 0.1 Ep to ^ Ep. The existence of non-zero su- 
perfiuid pair potentials induced by the SOC with a large 
Zeeman field can be detected through the emergence of 
vortices Q in experiments. The topological phase tran- 
sition can be observed by detecting the change of Eg 
versus F for each fixed kz using the momentum resolved 
photoemission spectroscopy 17, From Fig. UK, we 
see for l/Kpa, = -0.1 and F = 0.85^;^-, Eg ~ Ep, 
which is much larger than the energy resolution in the 
photoemission experiment. In experiments, the observed 
gap Eg first decreases to zero and then increases, which 
provides a clear signature for the transition from non- 
topological to topological superfiuid phases. Note that 
the transition temperature Tc from superfiuid to normal 
states in the presence of SOC should be larger than that 
for the superfluid without SOC [3| because of the en- 
hanced superfiuid order parameter A (see Fig. [2^) [lH. 
The temperature needed for the observation of the topo- 
logical phase crossover in Fig. 4b is clearly lower than Tc, 
which is nonzero on both sides of the critical Fc because 
of the existence of the non- vanishing A. 

Summary: In summary, we study the BCS-BEC 
crossover and the topological phase transitions in 3D uni- 
form spin-orbit coupled degenerate Fermi gases. The pre- 
dicted SOC induced s-wave superfluid opens new possi- 
bilities for generating and observing many new topologi- 
cal phenomena in Fermi gases. The observation of the 
3D topological phase transitions using the experimen- 
tally already realized photoemission spectroscopy pro- 
vides a critical first step for searching for non-trivial 
topological superfluid states (which in 2D support Majo- 
rana fermions and the associated non-Abelian statistics 
14, in cold atom s-wave superfluids, which are of 
not only fundamental but also technological importance. 

Note: After we submitted the paper, two manuscripts 
appeared in arXiv [2ll,[2^, where the BCS-BEC crossover 
with SOC (but without Zeeman fields) is discussed. 
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